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Abstract 

We apply the factorization technique developed by Kuru et. al. [Ann. Phys. 323 (2008) 413] 
to study complex classical systems. As an illustration we apply the technique to study the 
classical analogue of the exactly solvable VT symmetric Scarf II model, which exhibits the 
interesting phenomenon of spontaneous breakdown of VT symmetry at some critical point. 
As the parameters are tuned such that energy switches from real to complex conjugate pairs, 
the corresponding classical trajectories display a distinct characteristic feature — the closed 
orbits become open ones. 
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1 Introduction 



Fairly recently Kuru et. al. obtained exact solutions of the classical analogues of some exactly 
solvable one-dimensional quantum mechanical models in the framework of the factorization 
technique [1]. However, it must be kept in mind that this method cannot be applied to all 
classical systems as not all quantum mechanical models have exactly solvable classical analogues. 
On the other hand there have been several attempts to see if the spontaneous breakdown of VT 
symmetry at an exceptional point in a certain class of non Hermitian VT symmetric quantum 
mechanical systems gets manifested in the corresponding classical picture [21 HI O [6l El [8] . For 
example, the VT symmetric model H = p^ +x'^{ixY, was studied numerically in ref. [2] to show 
that the exceptional point occurs at e = in both the quantum as well as the corresponding 
classical systems. In a different work [6j the classical motion of two 1-dimensional non Hermitian 
systems was studied — viz., the complex Harmonic oscillator and the complex cubic potential 
— and an attempt was made to find a connection between the reality of the spectrum and 
the regularity of the classical trajectories. Motivated by these studies, in this work we intend 
to extend the factorization technique of ref. [1] to study complex classical systems, which 
are classical analogues of non Hermitian quantum mechanical systems. In particular, we shall 
apply the technique to study the classical analogue of the following exactly solvable quantum 
mechanical VT symmetric (complex) Scarf II potential [9j 

V{x) = —vi sech^x — -Oi sech x tanhx ; vi real , vi > (1) 

For vi pure imaginary (say vi = iv2), the model belongs to the well known category of VT 
symmetric potentials of the form V{x) = viVji{x) + iv2Vi{x) , with Vr{x) even, and Vi{x) 
odd. It is one of the few exactly solvable quantum mechanical models which exhibits the in- 
teresting phenomenon of spontaneous breakdown of VT symmetry at a critical value of the 
coupling parameter V2. Additionally, its classical analogue can be solved exactly by means of 
the interesting factorization method developed in ref. [T]. Our aim here is to see if the above 
system contains periodic or irregular trajectories, and check if the spontaneous breakdown of 
VT symmetry switching energy from real values to complex conjugate pairs, manifests itself in 
the corresponding classical system as well, in terms of non-periodicity of orbits or any other 
feature. 

The organization of the paper is as follows. To make the paper self contained, in Section 2 we 
briefly discuss the factorization technique of ref.[l] to obtain the classical trajectories of some 
exactly solvable one dimensional systems, and extend the same to complex classical systems. 
Based on the equations obtained in Section 2, in Section 3 we obtain the equations of the 
classical trajectories, the classical momenta and the phase-space trajectories of a particle under 
the influence of the complex Scarf II potential, and plot the corresponding figures for different 
values of the coupling parameter vi. Finally, Section 4 is kept for conclusions and discussions. 

2 Formalism 

We start with the one dimensional classical Hamiltonian (in units h = 2m = 1) 

H{x,p) =p'^ + V{x) (2) 

where x,p are the canonical coordinates, V{x) is the potential, and the Poisson bracket of x 
and p is given as {x,p} = 1. The equations of motion of the classical particle are given by the 
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Hamilton's equations 



dH 
dp 



2p 



P 



dH 

dx 



-V'{x) 



(3) 



Thus X = 2p = 2V'{x) which on integration gives the velocity of the particle as 

dx 



dt 



±2JE-V{x) 



(4) 



E being the energy of the classical particle. In our discussions since the particle is moving under 
the influence of complex forces, time t is treated as a real variable, and the path x{t) traced 
out by the particle as well as its velocity v can take complex values. The roots of the equation 
E — V{x) = give the locations of the turning points, while the initial conditions determine the 
initial velocity of the particle. It may be mentioned that by conventional notions the velocity of 
the particle can take real values only. So the initial positions of the particle under a real potential 
lie on the real axis, in between the turning points. However, if we analytically continue into the 
complex plane, a particle under the influence of a complex force can move about in the complex 
plane. Thus any point in the complex plane may be an initial starting point for the particle. 

Our aim here is to obtain exact analytical solutions of the orbits traced out by the classical 
particle. For this we shall mainly follow the factorization technique of ref. [Tj. We assume a 
factorization of the Hamiltonian H in the form 



H = A+A- + 7(iJ) 



(5) 



where unlike in usual quantum mechanical factorizations, ^{H) may depend on H, and A^ (not 
necessarily complex conjugate) are taken to be of the form 

A^ = ^if{x)p + VHg{x) + ^{x) + (l){H) (6) 

Furthermore, the functions A^ and H are assumed to define a deformed algebra with Poisson 
Brackets as follows : 

[A'^^H] = ±ia{H)A^ 

(7) 



where the auxiliary functions a{H), f3{H) and (j){H) are expressed in terms of the powers of 

In case the quantum version admits bound states with negative energies, ^/H should be replaced 

by \/—H. Making use of equations ([3]), 1^ and ([7]) we arrive at the following expressions : 

9 

fix) 



a{H) 



ip'{x) + g'{x)V H 



f{x)V'{x) - 2f'{x) [H - V{x)] = a{H) [g{x)VH + ^{x) + (t>{H) 



2Vh [nx)gix)-fix)g'ix)] 



g{x) [2nx)V{x) + f{x)V\x)] 



+ ^f{xf-^[H -V{x)]-2f{x) 



^ [x) + V (x) 



(8) 
(9) 

(10) 
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Now we construct two quantities of the form 



Q± = ^±e^^"(^)* (11) 
which are time dependent integrals of motion. Nevertheless, their total time derivative vanishes 

Thus 

I Q+Q- I = I A+A- I (13) 

so that the particular values of the integrals of motion and those of A^ may be respectively 
denoted by 

Q± = c(^)e^*^« (14) 
= c(S)e^*^^«+"(^)*^ (15) 
where is determined from initial conditions, and 



c{E) = ^E-^{H) (16) 

For c{E) to be real, the expression within the square root sign must be positive. This condition 
gives the range of energy values for the classical particle. The solution of (jl4p gives the trajec- 
tories x{t) and momenta p{t) of the corresponding classical particle in the complex plane. We 
illustrate our formalism with the help of an explicit example in the next section. 

3 Classical analogue of Complex Scarf II potential 

The quantum version of the "PT-symmetric (Complex) Scarf II potential displays certain inter- 
esting features — 

(i) its discrete spectrum is real below the VT threshold, often referred to as the Exceptional 
(or critical) point, whereas above it the energy values occur as complex conjugate pairs 

(ii) its continuous spectrum displays spectral singularity at the critical point, where the reflec- 
tion and transmission coefficients tend to diverge. 

For the sake of comparison of classical trajectories with other numerical studies (below and 
above the critical point), in this work we shall restrict ourselves to bound states only, hence, 
negative energies (i? < 0). The final expression for V{x) should be of the form 

V{x) = —70 sech^ + 26 sech tanh (17) 

The parameter 5 plays a crucial role here — the potential in (|17|) is real for real values of (5, VT 
symmetric for pure imaginary values of 5, and a general complex potential (without any VT 
symmetry) for complex values of 6. Our aim in this section is to find the exact equations for 
the classical trajectories and momenta of a particle moving in the complex plane under such a 
potential for different 5, with special emphasis on pure imaginary values, and to see if there is 
any change in these classical trajectories and momenta at the onset of spontaneous breakdown 
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of VT symmetry. The form of the potential in eq. (jl7p above demands that in the expression 
for in ([6]) we take the following forms of the functions g{x) , ^p{x) , 4>{H) : 

g{x)^0 , 99(x)=0 , ^(H) 



so that reduce to 

= Tif{x)p + V^g{x) + -j= (18) 
Solving equations ([8]), Q and (fTOj) simultaneously gives 

5f(x)=sinh— , /(x)=cosh— , 'y{H) = -jq + — (19) 



So using equations and p6j) . and the expression 
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A ' ^ = + 70 - 

we obtain the value of c{E) as 



A+A- =H + jo-- (20) 



c(i?) = y^ + 7o-^ (21) 

Equation ()2ip gives the range of values for energy E as c{E) should be real. It is worth remem- 
bering here that we are considering the energy E to be negative. We investigate the different 
cases in some detail below : 



Case 1 : (5 is real; 6 = 5r (say) 

This gives the real Scarf II potential, dealt with in ref. pLj, with classically allowed values of 
energy E lying in the range 



'"°'V"° +'^« < E < (22) 



Case 2 : 5 is pure imaginary; S = iSj (say) 

This is the classical analogue of the famous non Hermitian yet VT symmetric Scarf II potential, 
which exhibits the phenomenon of spontaneous VT symmetry breaking at the phase transition 
point, and real energy values switch to complex conjugate pairs. Straightforward algebra shows 
that for the corresponding classical motion c{E) reduces to 




c{E) = \E + ^o + -^ (23) 



Since we are dealing with negative energies, the classically allowed range for E in this case 
becomes 

-70- V^7o -4.5? -70 + V^7o -45/ ,,,, 
-7, < -E < (24) 
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One can check that the classical Hamiltonian, though complex, is symmetric under parity- 
time reversal. Furthermore, for 70 > | 26 j \, E is real, corresponding to exact or unbroken 
VT symmetry. However, for 70 < | 2(5/ |, energies turn complex. Thus the classical system 
undergoes a phase transition at 70 = | 25i \ . The corresponding quantum version, too, undergoes 
an abrupt phase transition from exact VT to spontaneously broken VT phase at some critical 
value of the coupling parameter 61. Our aim in this work is to see whether this phase transition 
gets manifested in the trajectories traced out by the classical particle. 

Case 3 : (5 is complex; 5 = 6ji + i5i (say) 

It is easy to observe that the potential (I17p is no longer VT symmetric; nor is it r/-pseudo 
Hermitian, and the energy spectrum is, in general, complex. Naturally, it does not arouse our 
interest; so we shall not pursue this case any further. 



3.1 Classical Trajectories 



To plot the classical orbits, we need to derive the expressions for the trajectories of the classical 
particle in the complex plane, under the influence of the potential (I17p . For this purpose we 
assume the most general form of c{E), viz., c{E) = cr{E) + icj{E) 

Putting 6 = 6r + i6i in the expression for , i.e., = c(£;)e=^*(''o+"o^^*), we obtain 



^ip cosh — — h 



— . aox 6r + i6i 
-h smh — — I — 



2 2 ■ y/-E 

= cr{E) cos (j)o + aoV-Et^ =F ci{E) sin (^60 + a^^f-Et 

+ i \ci{E) cos (6*0 + Qo\/^t) ± cr{E) sin (6*0 + ao\/^i) } 
and consequently, the final form of the trajectory x{t) and momenta p{t) as 

2 . \5r-cr{E)^J^ cos {eQ + aQ^/^t)±CI{E)^/^ sin {9Q + aQ^/^t 



(25) 



x{t) = — sinh 



E 



(26) 



E \^—cr{E) sin 

+ {Sr - cr{E)^/^ cos (6*0 + ao\/^t) ± CI{E)^/^ sin (Oq + ao\/^t) } 

(27) 

Equation (j26p describes the position and equation (j27p the momentum of a classical particle in 
the complex plane, as it is under the influence of complex forces. x{t) and p{t) given above are 
complex, and therefore we shall plot the real and imaginary parts of the position and momentum 
separately, to gain an insight into the behaviour of the particle's motion. In this work we shall 
primarily concentrate on purely imaginary values ol 5 = iSj, to see if any connection exists 
between the reality of the spectrum and regularity of the classical orbits. 
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3.2 VT symmetric Scarf II potential : 5 = i 5i 

Initially we start with bound states with real negative energies, in the unbroken VT phase in 
the quantum version. The corresponding classical condition for real energies is 2 | 5/ |< 79. To 
plot the exact trajectories it is essential to have a knowledge about the initial condition x(0). 
The classical turning points are given by the roots of the equation E — V{x) = 0, for which 
we need the precise value of E. The turning points are obtained in the form x± = ita + 
and hence are symmetric with respect to the imaginary axis. This is expected as the classical 
turning points of a complex T^T-symmetric potential for real energies, E, essentially occur 
as {z,—z*) : (—a + i6, a + ib) [9]. We plot the orbits traced out by the classical particle in 
the complex plane, in Fig. 1, for parameter values uq = 2 ,6 = 2i , 70 = 6. Equation 
(p^ suggests that classical energies should lie between —0.763932 and —5.23607. We plot the 
trajectories for E = —3. Any other value of energy within the specified range gives similar results. 
The classical turning points for this particular set of parameters are {±0.781368 — 0.528945z}, 
{±0.781368 - 2.61265i}, etc., i.e. of the type {zi,-z{) , (z2,-4) , etc. 

Closed Orbits for PT symmetric Scarf II model for unbroken PT symmetry 




-1.0 -0.5 0,0 0,? 1,0 



Re x(t) 

Figure 1: The classical trajectories for parameter values i? = — 3, ao = 2, (5 = 2z,7o=:6 
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In case the classical particle begins its motion on the real line, within the turning points or 
from any one of the turning points itself, say zi, the particle can only oscillate between the 
turning points {zi, —zl). Similar is the picture for the particle starting from any position on the 
real line within the turning points {z2, —Z2), or from any one of them. These paths are shown 
by the horizontal lines in Fig. 1. A trajectory joining any other pair of turning points (e.g. 
Zl, — ^2 or Zl, Z2 etc.) is forbidden because the particle is under the influence of a VT symmetric 
potential. However, if the particle starts its motion from any other point in the complex plane, 
the oscillatory trajectories are surrounded by closed orbits of a definite periodicity. As seen 
in Fig. 1, the trajectories are closed curves, which do not cross. The time period for each of 
the orbits in Fig. 1 is calculated to be 3.6276. Thus our observations for this exactly solvable 
analytical model are similar to those obtained in ref. [?,12J by numerical methods. In Fig. 2, 
we plot the momenta of the classical particle in the complex plane, and in Fig. 3 we plot the 
phase space curves, for the same set of parameter values as in Fig. 1. In these cases, depending 
on the initial starting point, the momenta plots as well as the phase space plots are closed 
curves without any crossings. However, while the trajectory plots and the phase space curves 
are symmetrical with respect to the imaginary axis, the momenta curves are symmetric with 
respect to the real axis. 



Momentum Plots for PT symmetiic Scarf II model for unbroken PT symmetry 




-0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Re p(l) 



Figure 2: The classical momenta for parameter values E = —3 , ao = 2 , S = 2i ,jo = 6 



Now if the parameter values are such that 2 \ 61 \ > 70 , then equation (f24l) suggests 
that the energies turn out to be complex conjugate pairs, in spite of the Hamiltonian being 
VT symmetric. This is referred to as the spontaneously broken VT symmetric phase in the 
quantum version. The classical trajectories are plotted in Fig. 4, for parameter values ao = 
2 , 5 = 2i ,7o = 3. From eq. (jl6p . the energies for this set of parameter values should lie 
between -1.5 + 1.32288i and -1.5 - 1.322881 We have taken E = -1.5 - 0.3i for the plots 
in Fig. 4. The abrupt phase transition in the quantum version gets manifested in the classical 
motion, too. Irrespective of the starting point of motion, all of a sudden the closed orbits 
become open and the trajectory loses its periodicity. The turning points are no longer of the 
type {z,—z*), and symmetry is lost. For this particular set of parameters the turning points 
occur at (-0.102199-0.470998^ , 0.102199- 2.67059i , -1.40526- 1.3489i , 1.40526- 1.79269i) 
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Phase Space Orbits for PT symmetric Scarf II model for unbroken PT symmetry 




-0.5 0.0 0.5 

Re x(t) 



Figure 3: The real part of phase space trajectories for parameter values i? = — 3 , — 2 , J = 2i , 70 = 6 
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etc, i.e. are of the form (a — ibi , —a — 162 ,c — idi 1 —c — id2)- The momenta plots and phase 
space trajectories, too, lose their regular behaviour — the curves are no longer closed ones, nor 
do they show any definite periodicity. 

4 Conclusions 

To conclude, we have extended the factorization technique of ref. [T] to study complex classical 
systems. As an illustration we have applied the technique to study the exactly solvable classical 
analogue of the exactly solvable, non Hermitian yet VT symmetric Scarf II model 

V{x) = —vi sech^x — i V2 sech x tanhx ; wi > 0, vi real 

Since the particle is under the influence of complex forces, it is expected to move about in the 
complex X plane. The quantum version of this particular Hamiltonian exhibits an abrupt phase 
transition at some critical value of f2, switching energy values from real to complex conjugate 
pairs. On the other hand, for complex values of ^2, there is no space-time reflection symmetry, 
and energies are always complex. We have given special emphasis on real values of V2 in this 
work, so as to study the effect of spontaneous VT symmetry breaking on the classical trajectories 
and momenta. Employing the factorization technique of ref. [T], we have found the equation of 
the path the classical particle traces out in the complex x plane, and plotted the corresponding 
orbits in Fig. 1 and Fig. 4. It is observed that so long as VT symmetry is unbroken, the 
motion of the corresponding classical particle is bounded, with all trajectories having the same 
periodicity — Fig. 1. Furthermore, the orbits are symmetric with respect to the imaginary 
axis, and the classical turning points are of the form (z, —z*). The spontaneous breakdown of 
VT symmetry resulting in switching of energy values from real to complex conjugate pairs, has 
an interesting manifestation in the corresponding classical picture — the closed periodic orbits 
abruptly become irregular and open, as shown in Fig. 4. However, none of the trajectories cross 
each other, as is expected for systems symmetric under VT- 

We have also plotted the momenta of the classical particle in the complex plane in Fig. 2, 
and the phase space trajectories in Fig. 3, so long as VT symmetry is unbroken and energies 
are real. Here, too, the curves are closed, and of definite periodicity. While the trajectory plots 
and phase space plots are symmetric with respect to the imaginary axis, the momenta plots 
are symmetric with respect to the real axis. At the onset of spontaneous breakdown of VT 
symmetry at the exceptional point, each of these plots abruptly loses its regular pattern, and 
the closed curves become open. 

To summarize, the observed deviation from regularity in the classical orbits at an exceptional 
point, for the analytical model considered in this work, is similar to that shown in earlier works 
O |6]. However, we follow the factorization approach given in [T] and thus our technique is 
different from the numerical / perturbative studies done in ref. [21 E]. Additionally, the model 
considered here is exactly solvable, both in its quantum as well as classical versions. 
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open orbits for PT symmetric Scarf II model for spontaneously broken PT symmetry 

0.5 I < < < < ^ < < < < ^ < < < < ^ < < < 



0.0 - 



-0.5 - 



-1.0 - 




E 

-1.5 - 



-2.0 - 



-2.5 - 



-3.0 - 




-0.5 0.0 0.5 

Re x(ll 



Figure 4: The classical trajectories for parameter values E = —1.5 — .3z , ckq = 2 , (5 = 2i ,70 = 3 
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